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After recapitulating the covariant formalism of equilibrium statistical mechanics in special rel- 
ativity and extending it to the case of a non-vanishing spin tensor, we show that the relativistic 
stress-energy tensor at thermodynamical equilibrium can be obtained from a functional derivative 
of the partition function with respect to the inverse temperature four-vector /3. For usual ther- 
modynamical equilibrium, the stress-energy tensor turns out to be the derivative of the relativistic 
thermodynamic potential current with respect to the four-vector /3, i.e. T M " = — <9<I? M /c>/3„. This for- 
mula establishes a relation between stress-energy tensor and entropy current at equilibrium possibly 
extendable to non-equilibrium hydrodynamics. 



One of the best known formulae of statistical mechan- 
ics states that, in the grand-canonical ensemble, the mean 
value of energy can be obtained as a derivative of the 
logarithm of the partition function with respect to the 
inverse temperature. In terms of densities, the energy 
density is the derivative of the pressure with respect to 
the inverse temperature, i.e.: 
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H/T 



In relativity, the mean value of energy density is the time- 
time component of the stress-energy tensor, so the ques- 
tion arises whether it is possible to obtain the full stress- 
energy tensor, and not just its time-time component, as 
a sort of derivative of the "density" of partition function. 
In this work we provide an answer to this question and 
we will show that: 



(1) 



that is the stress-energy tensor can be obtained as a vari- 
ational Euler-Lagrange derivative (in a sense which will 
become clear later) at thermodynamical equilibrium of 
the relativistic generalization of the density of the ther- 
modynamic potential log Z, a vector current $ M which is 
linearly related to the entropy current. 

In other words, equation ([T]) establishes a relation be- 
tween the mean value of the stress-energy tensor and the 
relativistic entropy current at equilibrium. If this rela- 
tion could be extended to non-equilibrium situations, one 
would have a tool to determine entropy current from the 
expression of the stress-energy tensor, which is one of the 
main problems of dissipative relativistic hydrodynamics 
0, [1] (for a recent discussion see also ref. Q ) ■ The ex- 
pression of entropy current in non-equilibrium is actually 
used as means to determine the structure of the stress- 
energy tensor itself and the generalization of the relation 
([1]) would therefore be very important. 

Before proving eq. ([T]), we will recapitulate the fully 
covariant formulation of equilibrium in relativistic sta- 
tistical mechanics, including the possibility of a non- 
vanishing spin tensor, hence of a non-symmetric stress- 
energy tensor. We start with a brief summary of equi- 



librium thermodynamics of a systems with angular mo- 
mentum, which macroscopically corresponds to a rigidly 
rotating fluid [4|]. A detailed discussion of this system 
can be found in refs. [5|. 

Thermodynamical equilibrium occurs when entropy is 
maximal. The maximization of S = — tr(plogp) with 
the constraint of fixed, constant mean energy, charge and 
mean total angular momentum leads, as it is well known, 
to the density operator 0, Q : 
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P : 



exp[-H/T + luJ z /T + fj,Q/T] 
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where Z = tr(exp[-ff/T + ojJ z /T + fiQ/T]). The phys- 
ical meaning of T is that of temperature of an energy 
reservoir and uj that of an angular velocity of an angular 
momentum reservoir in contact with the system. In a 
quantum relativistic field theory, the operators in eq. ([2]) 
can be written as integrals over some space region V: 
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S being the spin tensor. The operator ([2]) can be written 
in a fully covariant form. First define: 



= (i/r)(i,wxx) 
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with u) = wk constant angular velocity vector directed 
along the z axis. Then: 



uj Xv = w/T(SlS 2 v -5iSt) 
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which is the acceleration tensor for a rigid rotation 
Q. Finally, define the normal versor of the three- 
dimensional space-like hypersurface V (embedded in 
Minkowski spacetime) appearing in the eq. as t and 
its measure dS M = d 3 xt^. Hence, we can rewrite ((2|) as: 



9 = ~z exp 



dS, 



The four-vector /3 in eqs. ((4]) is then, by construc- 
tion, the inverse temperature four- vector and 1 / \[W = 
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T I y/l — \\u> x x|| 2 = To is the invariant temperature, 
i.e. the temperature measured by a thermometer mov- 
ing with the rigid velocity field u> x x with respect to 
the thermostat at temperature T [JJ, Q ■ The latter ex- 
pression does not fulfill yet the request of full covariance 
as it apparently depends on a particular hypersurface V. 
In fact, if the divergence of the integrand vanishes and 
if its flux at the boundary of V, i.e. dV, also vanishes, 
the space-like hypersurface is arbitrary and the density 
operator can be finally written as: 
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where E is a general, arbitrary, space-like hypersurface 
bounded by the same dV. The covariant form ^ of 
the equilibrium statistical operator, to our knowledge, 
was first written down by Weldon Q; the above form 
generalizes his formula in that it includes a non- vanishing 
spin tensor, which is generally needed if the stress-energy 
tensor is not the Belinfante symmetrized stress-energy 
tensor ||. 

The two aforementioned conditions on the integrand 
also ensure the stationarity of the density operator with 
respect to any inertial frame, because p is unchanged in 
© going from the hypersurface V(t) to V(t + At). The 
divergence of the integrand reads: 



9a -T^f3„ 
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where we have taken into account the continuity equa- 
tions of the stress-energy tensor and the current. As 
£ = n/T and w\ v (see eq. [5} are constant, the above ex- 
pression reduces to, after separating the symmetric and 
antisymmetric part of the stress-energy tensor: 
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Because of the continuity equation of the angular mo- 
mentum tensor: 



the 



vanishes if: 



d u & M + 2T\ V = 



d»Pu + 9„/3 u = 



(9) 
(10) 



It can be readily checked that the four- vector (3 and the 
tensor w in eqs. (jj) and §5§ fulfill both conditions above. 
Hence, the divergence of the integrand does vanish: 



9 U [ -T^ft, + Is^uxu + ft" I = 
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The second condition, i.e the vanishing at the boundary: 

1 
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(-T + ^ Xv u Xv + (12) 

must be enforced through suitable boundary conditions 
of the quantum fields. 

Instead of the equilibrium density operator ([2]) , we can 
use the general covariant formula ((6]) as a starting point 
and look for the implications on /?, w and £ of thermo- 
dynamical equilibrium conditions, i.e. of eqs. (|11[) and 
(|12[) ; in principle, with this approach, we could find new 
forms of equilibrium distributions. For the divergence to 
vanish, according to eq. ©, one needs to have: 



9xt = 



9 X u U v = 



and, again, the equations (|10l) . The first of (|10l) entails 
that the inverse temperature four-vector ought to be a 
Killing vector, a well known condition for equilibrium 
[H, 0. Together with the second of eqs. (ITU1) . this leads 



(13) 



with b constant four-vector. This expression of the 
inverse temperature four-vector comprises all possible 
forms of relativistic thermodynamical equilibria; the ro- 
tating case ((U) actually corresponds to 6 M = (1/T, 0) and 
u) given by eq. ([5]). 

By using the eq. ([6]) one can write down the entropy: 
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S = logZ+^d^ (t» v $ v - ~S»> Xv u; Xv - 

where the symbols without hat denote the mean values of 
quantum operators (i.e. tr(pA) = A) 0. For an entropy 
current to exist in relativistic thermodynamics [l|, [2J so 
that the total entropy can be written as S — J s dE M s M , 
the logarithm of the partition function must be written as 
an integral over the same hypersurface E of a vector field 
tfr^, hereby defined as thermodynamic potential current: 



logZ 



dE M ^ 



so that the entropy current reads: 

s " = $m + T^p v - -S^ Xv uj X u - ft 



(15) 
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1 It should be pointed out that the mean value of operators in- 
volving quantum relativistic fields are generally divergent (e.g. 
T 00 for a free field has an infinite zero point value). To remove 
the infinities, the mean values must be renormalized, what can 
be simply done for free fields by using normal ordering in all 
expressions, including the density operator itself. Henceforth, it 
will be understood that all the mean values of operators are the 
renormalized ones. 

2 Note that in this formula log Z is meant to be the renormalized 
one, see previous footnote 
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In principle, the existence of the thermodynamic poten- 
tial current could be proved working out the trace of 
Zp in eq. ([6]) , but this requires lenghty manipulations of 
the commutators of T^" and so we rather assume (|15j) 
invoking the existence of an entropy current. At equilib- 
rium, the four-vector field $ must be divergence-free and 
it must have vanishing outward flux through <9E for the 
partition function ()15[) to be stationary and independent 
of the space-like hypersurface E. In view of the p^|) , the 
condition 3^$^ = is indeed a consequence of the re- 
quirement of vanishing entropy production (d^s^ = 0) at 
equilibrium. We point out that the thermodynamic po- 
tential current is not uniquely defined, as one may add 
a divergence of an antisymmetric tensor field with suit- 
able boundary conditions to obtain the same partition 
function in eq. (|15|) . 

The mean value of the stress-energy tensor at equilib- 
rium reads: 



T^ix) = tr(pf^) = -tr ( T^(ir)exp 

z 



dE, 
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where © has been used. Let us now fix the space-like 
hypersurface E and write dE^ = dEn M , where is its 
normal time-like unit vector. One can obtain a contrac- 
tion of the mean value of the stress-energy tensor with 
the normal vector by taking a functional derivative 
with respect to the four-temperature vector j3 seen as a 
function of x, keeping oj and £ fixed, with respect to the 
measure dE: 
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Formally, the above formula can be shown by using the 
expansion of the exponential of sum of operators (Zassen- 
haus formula) and taking advantage of the ciclicity of the 
trace. While the left hand side depends on a vector n M , 
which is arbitrary, the rightmost hand side is not mani- 
festly dependent on it. In fact, the functional derivative 
of the partition function includes a hidden dependence 
on the normal vector as the functional derivation implies 
the choice of a measure, hence of a hypersurface E and a 
corresponding normal vector. 
In view of the formula (IT5|) : 
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At equilibrium, the thermodynamic potential current de- 
pends on the equilibrium values of /3,ui and £. Taking the 
functional derivative means moving the field j3 slightly 
out of equilibrium, i.e. P{x) = /3 eq (x) + 5/3(x), what may 
introduce dependences of <f> M on the derivatives of the /3 



field. As it is well known from the theory of functional 
derivation, provided that the perturbation 6/3 is chosen 
so as to fulfills suitable boundary conditions, the (Til?]) 
yields: 



rvT^(x) = n 
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where . . . stands for terms involving higher order deriva- 
tives of the P field. However, since the vector n^ is an 
arbitrary timelike field, the straightforward consequence 
of (HDJ is: 



rj-iflU 
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i.e. the stress-energy tensor is minus the Euler-Lagrange 
derivative, denoted by $ and defined by eq. (|20p , of the 
thermodynamic potential current. 
Similarly, it can be shown that: 



(22) 



cq 



We are now going to work out the above formulae in the 
simplest instance of thermodynamic equilibrium, which 
is the familiar one with — 1/Tqu^, with = const 
normalized four- velocity, £ = u/T = /Uq/Tq = const , 
and u> = 0, i.e. no rotation. Thus ^ = = l/T u^ 
with To the invariant temperature by definition. The 
resulting density operator p = 1/Zexp[—P ■ (3 + Q£\ is 
invariant for translations, implying that all mean values 
of fields (including stress-energy tensor) are constant in 
space-time. The thermodynamic potential current $^ at 
equilibrium is a vector function of fi, uj and £, but since 
ft is the only non-vanishing vector field, it can only be of 
the form: 



(23) 



where the physical meaning of the scalar function p(/3 2 , £) 
is to be found. The above form of the thermodynamic 
potential current is unambiguous because, due to the con- 
stancy of the arguments, any additional divergence of an 
antisymmetric tensor field vanishes. For this special kind 
of equilibrium, it was written down first in ref. Q ■ 

In eq. (|2"Tj) . all terms of the Euler-Lagrange derivative 
of the $ function involving derivatives vanish at equilib- 
rium because, e.g.: 
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and since all derivatives of the f3 field vanish at equilib- 
rium, this term altogether vanishes. The same applies to 
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possible terms involving derivative of the $ function with 
respect to higher order derivatives of the /3 field. Further- 
more, it is not difficult to realize that the only term of 
(9$ A1 /d(3 u contributing at equilibrium is the derivative of 
(|23| itself, that is: 



d/3 v 



d<S>v 
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and, similarly: 



eq 



dp 



/?" = nu>* 
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(26) 



where (|23|) has been used. By using ([25]) and ([231 we get: 



dp 
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(27) 



This form would be enough to identify p as the pressure, 
because if u = (1, 0) then p is the diagonal element of the 
spacial part of the stress-energy tensor. This identifica- 
tion is confirmed by the expression of the proper energy 
density p obtained from ([2"T]) : 



and the eq. (|27|) turns into the familiar: 

T" v = (p + p)u^u v - pg^ 



(28) 



(29) 



The equation ([28]), that is 2dp/d/3 2 \ i = ~(p + p)/(3 2 is 
just an alias of the Gibbs-Duhem relation dp/dT^^ = s 
which can be readily checked taking into account that 
the above derivative is taken by keeping £ = pq/Tq fixed 
and that Tqs = p + p — p^n (see below). Similarly, the 
eq. (|2"oT) is an alias of the relation dp/dp \T — n. 

Finally, we show that all known thermodynamic rela- 
tions involving proper entropy density are also recovered. 
From (p~6|) with ui — 0, using ([23]) we get, by contracting 
with the four- velocity u: 



v = v\W + p\/W - i n 



(30) 



which, by using — 1/T and £ = Po/Tq reads as 

the familiar relation Tqs = p + p — p$n. We can obtain 
the differential of cntroopy current (difference between 
nearby equilibrium states with /3 and £ being the under- 
stood parameters): 

d s ^ = + T^&fiv + /3„dT^ - £dj>* - jM£ (31) 



Since: 



in view of the J25 



26J), the eq. 



d£ 



311) turns into: 



ds" = ^dT^ - £df 



(32) 



(33) 



This equation was obtained by Israel and Stewart re- 
sorting to several assumptions concerning the familiar 
form of thermodynamical equilibrium, including the ex- 
pression (|29p itself; in fact, in our method, both ([55]) and 
([29]) are consequences of the eq. (|2T|) . derived in turn 
from the general form of statistical density operator ([B]), 
which is more economical and transparent. Contracting 
the above formula with u, t we get: 



u^ds 11 =ds- s^du^ = ds - + T^P V - tf)du^ 



(34) 



The expression within round brackets is parallel to 
and this makes the whole term vanishing as u^du^ — 
(l/2)dw 2 = 0. For the same reason, one has u^u^dT^ = 
diu^UvdT^) = dp and £,u^df = £dn. Therefore, the 
eq. (|33[) becomes the well known: 



Tods = dp — podn 



(35) 



We stress that all formulae (|23ti35|) apply to the ther- 
modynamical equilibrium without rotation, i.e. with 
bj = in the @. The application of (|2"Tj) to the most 
general form of equilibrium (1131) . and to slightly out-of- 
equilibrium situations may shed light on a general re- 
lation between entropy current and the mean values of 
stress-energy, spin tensors and the charge current vector. 
This will be the subject of further studies. 
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